#Exerccio5 - Liedson Tavares de Sousa Carneiro

##5.1

esaligna<-read.table("esaligna.csv",head=T,sep=",",as.is=T)
str(esaligna)
head(esaligna)
#1
plot(esaligna$ht~esaligna$dap, xlab="Altura(m)", ylab="DAP(cm)")

#2 #3 #4
plot(esaligna$ht~esaligna$dap, xlab="Altura(m)", ylab="DAP(cm)", main="E. saligna", font.main=3, bty="l",tcl=0.3)

#5
par(cex=1.5)

##5.2
factor(esaligna$talhao)
boxplott<- data.frame(esaligna$talhao,esaligna$dap)
colnames(boxplott)<- c("talhao","DAP")
x11()
par(mar=c(14,4,8,2))
par(mfrow=c(1,2))
boxplot(boxplott$DAP~boxplott$talhao, xlab="Talhes", ylab="DAP(cm)")
text(x=0.6,y=23,"A")
meant<-tapply(X=esaligna$dap,INDEX=esaligna$talhao, FUN=mean)
meant
sdt<-tapply(X=esaligna$dap, INDEX=esaligna$talhao, FUN=sd)
sdt
talhao<-data.frame(meant,sdt)
class(talhao)
talhao2<-matrix(1:18,nrow=6,ncol=3)

plot(talhao$meant,axes=F,xlab="Talhes",ylab="Mdia DAP(cm)",ylim=(c(0,25)))
axis(1,at=c(1,2,3,4,5,6),labels=c("16","17","18","22","23","32")) 
axis(2,at=c(0,5,10,15,20,25))
par(new=TRUE)

plot((c(talhao$meant+talhao$sdt)), axes=FALSE, ann=FALSE, pch="-",ylim=(c(0,25)))
par(new=TRUE)
plot((c(talhao$meant-talhao$sdt)), axes=FALSE, ann=FALSE, pch="-",ylim=(c(0,25))) 
arrows(x0=c(1,2,3,4,5,6),y0=(talhao$meant-talhao$sdt),x1=c(1,2,3,4,5,6),y1=(c(talhao$meant+talhao$sdt)), code=0)
text(x=1.1,y=24,"B")

#5.3
xy<- read.table("exercicio3-5.3.csv",head=T,sep=",",as.is=T)
str(xy)
factor(xy$y2)
length(xy$x1)
x11()
par(mar=c(14,4,8,2))
par(mfrow=c(1,2))
plot(xy$y1~xy$x1,xlab="Log(Patch size)(ha)",pch=17,tcl=0.3,bty="l",font=6,ylab="Euclidean distances",xlim=(c(0.5,2.3)),ylim=(c(0.0,3.0)))
segments(0.5, 2.3 , 2.15, 0.97, lwd=2)
text(x=2.2,y=3.0,"a",cex=1.5,font=6)
boxplot(x2~y2,data=xy,axes=F,tcl=0.3,bty="l",xlim=(c(0,7)),ylim=(c(0.0,3.0)),xaxt="n",yaxt="n",outline=F)
axis(1,at=c(0,1,2,3,4,5,6,7),labels=c("","Small","Medium\nEdge","Medium\nInterior","Large\nEdge","Large\nInterior","Control",""),tcl=0.3,font=6)
axis(2,c(0.0,0.5,1.0,1.5,2.0,2.5,3.0),tcl=0.3,font=6)
par(new=TRUE)
text(x=(c(1,2,2.9,3.1,4,4.8,5,5.2)),y=(c(3,3,3,3,3,3,3,3)),"*",cex=1.5)
text(x=6.5,y=3,"b",font=6, cex=1.5)